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ABSTRACT 

Rotation curves of galaxies show a wide range of shapes, which can be paramaterized as scat¬ 
ter in Vrot(lkpc)A^rnax i-^. the ratio of the rotation velocity measured at 1 kpc and the maxi¬ 
mum measured rotation velocity. We examine whether the observed scatter can be accounted 
for by combining scatters in disc scale-lengths, the concentration-halo mass relation, and the 
M^-Mhaio relation. We use these scatters to create model galaxy populations; when housed 
within dark matter halos that have universal, NFW density profiles, the model does not match 
the lowest observed values of Vrot(lkpc)Wmax and has too little scatter in Vrot(lkpc)Wmax 
compared to observations. By contrast, a model using a mass dependent dark matter profile, 
where the inner slope is determined by the ratio of M^/Mhaio, produces galaxies with low val¬ 
ues of Vrot(lkpc)Wniax dud a much larger scatter, both in agreement with observation. We 
conclude that the large observed scatter in Vrot(lkpc)Wniax favours density profiles that are 
significantly affected by baryonic processes. Alternative dark matter core formation models 
such as SIDM may also account for the observed variation in rotation curve shapes, but these 
observations may provide important constraints in terms of core sizes, and whether they vary 
with halo mass and/or merger history. 
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1 INTRODUCTION 

The rotation curves of disc galaxies are diverse (e.g. |Zwaan et al.| 
|1995[[Swaters et al.|2009| ); even galaxies with similar velocities at 
the flat part of the rotation curve can have significantly different 
rotation velocities at small radii. This results in a large scatter in 
rotation velocity at small radii, when plotted as a function of Vmax 
( [Oman et al.|2015| who measure Vrot at 2 kpc). 

It is well documented that a population of galaxies have slowly 
rising rotation curves, meaning a low value of Vrot(lkpc) compared 
to Vmax. Such galaxics are presented as evidence for dark matter 
profiles that are relatively flat in the inner region, often referred to 
as cored profiles (|Moore|1994||Salucci et al.|2003||Kuzio de Naray| 
let al.|2008[|0h et al.|2011t . It is argued that this contradicts a pre¬ 

diction of cold dark matter cosmology, that dark mater halos have 
steep inner density profiles, or “cusps” ( [Navarro et al.|1996] ). Yet 
other galaxies have rotation curves that do rise steeply in the inner 
region, and which can be fit by cuspy dark matter density profiles, 
having relatively high values of Vrot(lkpc) compared to Vmax- It 
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would seem that a successful model of galaxy formation needs to 
account for both cored and cuspy density profiles, sometimes in 
galaxies that have similar dynamical masses. 


This Letter examines the scatter in the Vrot(lkpc)-Vmax plane 
in a model population of galaxies that incorporates the empirical 
scatter in disc scale-lengths ( jCourteau et al.||2007| jS waters et al.j 
|200^|Fathi et al.|2010|), scatter in the concentration-halo mass re¬ 
lation ( [Dutton & Maccio|2014| t, and scatter in the M*-Mhaio rela¬ 
tion ( [Reddick et al.[[2013| Behroozi et al.[[2013[ ). Different forms 
of the underlying dark matter halo are assumed. Firstly the univer¬ 
sal “NFW” density profile that provides a good fit to dark matter 
only simulations ( [Navarro et al.[[199'^ . Secondly a mass depen¬ 
dent density profile that accounts for the effects of baryons on dark 
matter ( [Di Cintio et al.|20T4^ ,b). And finally, a cored density pro¬ 
file that may occur in alternative dark matter theories, in particular 
self-interacting dark matter (SIDM). 


In Section 113 we present our model galaxy populations. In 
Section [2^ we compile an observational sample of galaxies which 
have detailed rotation curve data. In Section [3] we show the model 
populations in the Vrot(Ikpc)-Vmax plane, and plot the scatter in 
Vrot(lkpc)/Vmax as a function of Vmax, making comparisons with 
the observational sample. We discuss our findings in Section|^ con¬ 
cluding that baryonic physics play a role in shaping the dark matter 
density profiles in the inner 1 kpc of galaxies. 
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Figure 1. Scale length versus luminosity for 144 observed galaxies (Hunter 
& Elmegreen 2006), which is ht by the red line, around which there is a 
scatter of ~0.23 dex 


Figure 2. 18 rotation curves are shown for the models with NEW and DC 14 
density profiles in the left and right panels respectively. The total dark mat¬ 
ter halo masses associated with this subset of model galaxies have maxi¬ 
mum velocities 90 < V<100 km s “ ^. 


2 METHOD 
2.1 Model 


Our model galaxies consist of exponential stellar and gas discs em¬ 
bedded within dark matter halos. For the dark matter halos, we 
use the concentration-mass relation for Planck ACDM ( [Dutton &] 
|M accio l20T4t , including a correction for adiabatic contraction, as 
found in cosmological simulations ( |Di Cintio et al.|[2014| ). This 
contraction correction alters the concentration rather than the in¬ 
ner slope, allowing the inner slope of the density profile to be set 
by hand, and has the form: 

C = ( 1.0 + O.OOOOSe^'^^) X Cnfw ( 1 ) 


where X = log^o(^*/^haio) + 4.5, and Cnfw is the concentra¬ 
tion given from the adopted ( [Dutton & Maccio|2014| l relation. 

Galaxy stellar masses are matched to dark matter ha¬ 
los using the empirical abundance matching relation ( [Guo[ 
[et al.[[ 20 T 0 ] ). We include the universal dark matter fraction of 
mass (QDM/^lm=0.83) in dark matter halos, and add a stel¬ 
lar disc with scale-lengths (hg) taken from the observed M,^- 
hs relation, logio(hs)=—2.8-i-0.351ogio(M*), (Dutton 07 gives 
logio(hs)=—2.46-i-0.2811ogio(M*)) coming from a fit to galax¬ 
ies from [Hunter & Elmegreen[ ( [2006[ ), which is similar to 
the relation found previously in larger studies of generally 
more massive galaxies ( [Courteau et al.[ [2007[ [Fathi et al.[ 
[2010[ ). We also add a gas disk with mass Mgas=1.3MHi, with 
logio(MHi/M,t)=—0.431ogio(M*) + 3.75 ( [Papastergis et al. 20T^ 
and disc scale length h^=3hs. 

The mass dependent ‘DC 14’ density profile is based on cos¬ 
mological simulations which match a wide range of galaxy scaling 
relations ( [Brook et al.|2012[[Stinson et al.|2013[[Kannan et al.|2014[ 
[Obreja et al.[2014| ). The DC 14 profile accounts for the expansion of 
dark matter halos due to the effects of feedback from star forming 
regions ( [Pi Cintio et al.[2014] |, taking the form ( [Merritt et al.[2006[ ): 

Ps 


p(r) = 




«-| (/3-7)/o; 


( 2 ) 


where (a, /3, 7 ) are the sharpness of the transition, the outer and the 
inner slope, respectively. 

The DC 14 profile allows for a range of inner slopes, 
determined by the stellar-to-halo mass ratio. Galaxies with 
M*<3 X 10 ® M 0 are dark matter dominated, and do not produce 
enough energy to flatten the halo’s inner density profile, which re¬ 
mains steep ( [Penarrubia et al.|20T^[Governato et al.| 2012 |[Di Cin-[ 


[tip et al.[2014[ |. The inner density profile then becomes increasingly 
flat as stellar mass increases relative to the dark matter mass ( [Gov-[ 
[ernato et al.[2012[|Di Cintio et al.[2014[ l, with the greatest flattening 
when M,t~3x 1 O®M 0 . For higher masses, the deeper potential well 
is able to oppose the halo expansion ( [Di Cintio et al.|2014^ ,b), re¬ 
sulting in a profile which becomes steeper, and returns to the NFW 
value at about the Milky Way mass. 

We make a simple parameterisation as representative of the 
self interacting dark matter (SIDM) case, which is loosely based 
on [Zavala et al.[ ( [2013[ ). Core size is set to zero for Milky Way 
mass galaxies (M*= 6 x 1 O^®M 0 ) and increases as the scale-length 
decreases, until reaching a maximum for M* = 1 O®M 0 , with lower 
mass halos all having the maximum core size, for which we use 
values of 1, 2 and 3 kpc. Cored profiles for the SIDM model have 
the form: 


p{r) = 


PoTs 


( 3 ) 


(tc + r){rs + r)2 
where po is a characteristic halo density, Vg is a scale radius and Vc 
is the core size ( [Penarrubia et al.|2012[ ). 

Note that this definition of core size, set by the parameter Vc, 
does not necessarily match definitions used in studies which use 
different density profiles, so care needs to be taken when comparing 
core sizes quoted here to those quoted in other studies. 

The important distinction between the SIDM model and the 
DC 14 model is that in the SIDM case, core size depends on halo 
mass, while in the DC 14 model the core size depends on the ratio 
of M*/Mhaio- We have not included scatter in the relation between 
Mhaio and core-size which may occur in SIDM models, dependent 
on different merger histories. Such a scatter is apparent in high 
mass SIDM halos [Rocha et al.[ ( |2013| ), but has not been quantified 
for low mass SIDM halos. 


2.2 Scatter 

We then add scatter to our models. Scatter in disc scale length 
comes from the empirical scatter of 0.23 dex in the luminosty- 
disc scale length relation, determined using 144 observed galaxies 
( [Hunter & Elmegreen|2006[ Oh et al 2015 in prep), as shown in Fig- 
ure[T] This is lower than the scatter of ~0.3 dex found in [Courted 
[et al.[ ( [2007[ ) for a larger sample of late type galaxies, but a sam¬ 
ple which is of generally higher mass galaxies than used in this 
study. [Fathi et al.[ ( p010[ ) explored the relation between scale-length 
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Figure 3. Vrot(lkpc) versus Vmax for the model galaxy populations when 
assuming the NFW (top panel, blue diamonds) and DC 14 (bottom panel, 
red squares) profiles. Black circles are observed galaxies. Values for halos 
from dark matter only simulations, including Icr scatter, are shown as a grey 
band. The 18 galaxies shown in Figure|^are marked with ’+’ symbols. 


and stellar mass down to lO^'^M©, and found similar scatter to 
[Courteau et al.| ( [2007] l, but with lower scatter for low mass galaxies. 
We tested our results using a scatter of 0.3 dex and found relatively 
insignificant effects that do not alter our conclusions. Scatter in the 
mass-concentration relation is set at 0.11 dex ( [Dutton & Maccio 
2014 1. Finally, scatter in M*-Mhaio is set at 0.2 dex ( [Reddick et al. 
2013 [Behroozi et al.[2013] l. 


2.3 Observations 


We compiled 110 observed galaxi es ([de Blok & Bosm^[2002[ 


Swaters et al.[[2003| 

de Blok et al.[[2004[ 

Kuzio de Naray et al.[ 

20081 [de Blok et al. 

[20081 [S waters et al.[ 

20091 |Oh et al.||20111 


Mannheim & O’Brien[2012[ Oh et al 2015, and references therein) 
with Vmax<150kms“^ that have rotation curves that allow us to 
compare Vrot(lkpc) with Vmax- 


3 RESULTS 

In Figure we show 18 rotation curves for the models with NFW 
and DC 14 density profiles in the left and right panels respectively. 
This subset of model galaxies have associated dark matter halos 
with maximum velocities 90<Vj^^<100kms“^: the only differ¬ 
ence between the two panels are the density profiles. A signifi- 
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Figure 4. Standard deviation of the ratio Vrot(lkpc)Wmax versus Vmax 
for the NFW (blue diamonds) and DC 14 (red squares) models. Observed 
galaxies are shown as black circles. Error bars show the standard errors. 


cant range of rotation curve shapes is apparent in both models. The 
DC 14 models have rotation curves that are more slowly rising than 
the NFW models. 

To measure the diversity in rotation curve shapes we plot, in 
Figurej^ Vrot(lkpc) as a function Vmax for the model galaxy pop¬ 
ulations when assuming the NFW (top panel, blue diamonds) and 
DC 14 (bottom panel, red squares) profiles. Vmax is defined as the 
rotation velocity as measured at 5 disc scale lengths, better repre¬ 
sentative of the observational data than simply plotting maximum 
circular velocity. Also shown in each panel, as black circles, are the 
observed galaxies. The value of Vmax for halos from dark matter 
only simulations, including Icr scatter, are shown in each panel as a 
grey band. The 18 model galaxies whose rotation curves are shown 
in Figureare marked with ’+’ symbols. 

Two features of Figure are worthy of attention. 
Firstly, around half the observed galaxies have values of 
Vrot(lkpc)/Vmaxthat fall below the NFW model population. This 
is just re-stating that a significant number of observed galax¬ 
ies are better fit by cored profiles than cuspy profiles. By con¬ 
trast, the DC 14 model population includes galaxies with low 
Vrot(lkpc)/Vmax, better matching the observations. 

The second feature is that the models with universal NFW pro¬ 
files scatter around (and slightly above) the expected values from 
pure N-body simulations, with no apparent shape to the distribu¬ 
tion, in terms of mass dependance. This refiects the universality of 
the NFW profile. By contrast, both the observed galaxies and the 
DC 14 model galaxies appear to have a structure in their distribu¬ 
tion; a relatively flat relation at low velocities, with a positive re¬ 
lation between Vrot(lkpc) and Vmax for higher rotation velocities, 
i.e those with Vmax>70 km s“^. 

We further explore the scatter in Vrot(lkpc)/Vmax in Figure]^ 
where the standard deviation of the ratio Vrot(lkpc)/Vmax is plot¬ 
ted as a function of Vmax for the NFW (blue diamonds) and DC 14 
(red squares) models. The observed galaxies are again shown as 
black circles. In each case, error bars indicate the standard er¬ 
ror. Clearly the NFW model has far too little scatter in the ratio 
Vrot(lkpc)/Vmax at all values of Vmax- The DC 14 model is a bet¬ 
ter match to observations. The match is certainly not perfect, with 
the DC 14 model having a peak dispersion at a slightly lower value 
of Vmax than is found in the observed population. 
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Figure 5. The relation between and V(lkpc), plotted for the NFW (blue 
diamonds, top panel) and DC14 (red squares, bottom panel) models, with 
observations shown as black circles. 

3.1 M,,-V(lkpc) 

We now look at the relation between and V(lkpc), plotted 
for the NFW (top panel) and DC 14 (bottom panel) models in 
Fig. figimsvL It is clear from Fig. figimsvl that the NFtW model 
is not able to form galaxies with low values of V(lkpc) at all stel¬ 
lar masses. By contrast, the observed relation between stellar mass 
and V(lkpc) is well produced by the DC 14 model. We caution that 
the data set is not homogeneous in terms of stellar masses. A ho¬ 
mogeneous set of galaxy stellar masses, along with high resolution 
rotation curves, as being compiled by Lelli et al. 2015 (in prep), 
would be invaluable here. 


3.2 Self Interacting Dark Matter 

We made similar plots for the SIDM model population. Figure 
shows Vrot(lkpc) as a function of Vmax as green triangles, where 
we have used a maximum core size of 3 kpc, with observed galax¬ 
ies over-plotted as black circles. Again, expected values for dark 
matter only simulations, including la scatter, are shown as a grey 
band. The galaxies associated with dark matter halos with maxi¬ 
mum velocities 90<V™<100kms“^ are marked with ’-i-’ sym¬ 
bols. Figurej^shows the scatter in the ratio Vrot(lkpc)/Vmax plot¬ 
ted as a function of Vmax for the SIDM model for 3 different core 
sizes, 1 kpc(pink triangles), 2kpc(purple triangles) and 3 kpc(green 
triangles) compared to the observed galaxies. 



Figure 6. Vrot(lkpc) versus Vmax for the SIDM model population with 
core size 3kpc(green triangles) and observed galaxies (black circles). The 
expected values for halos from dark matter only simulations, including la 
scatter, are shown as a grey band. Galaxies associated with dark matter ha¬ 
los with maximum velocities 90<V^^<100kms“^ are marked with ’+’ 
symbols 

For the paramaterization that we have used, the SIDM model 
requires a relatively large core size, 3 kpc, in order to reproduce the 
observed relations. As is evident in Figurej^ SIDM model galaxies 
can have low values of Vrot(lkpc), at given Vmax, as cored profiles 
make a better match for galaxies with slowly rising rotation curves. 
The SIDM models also have a larger range of rotation curve shapes 
than the NFW models. The larger the cores in SIDM halos, the 
greater is the effect of the scatter in the baryonic mass distribution 
(i.e. scatter in scale-lengths) in determining the variation in rotation 
curve shapes. The SIDM model with 3 kpccores is able to do as 
well as the DC 14 model in reproducing the observed scatter in the 
rotation curve shapes. 

What has not been included in the SIDM model is a depen¬ 
dence of core size on mass accretion history, which would result in 
scatter in the core size at given halo mass. This is expected, and has 
been shown on cluster scales |^cha et al.H2013| ), but has not been 
quantified at low masses. This scatter will likely make the SIDM 
model closer to observations, and difficult to distinguish from core 
formation through baryonic processes, although the detailed rela¬ 
tion between scatter and galaxy mass may still differ between the 
models. 


4 CONCLUSIONS 

We have created a model population of disc galaxies, using empir¬ 
ical relations between disc scale length and stellar mass, between 
stellar and gas mass, and between stellar and halo mass. Scatter in 
the M*-scale length and M^^-Mhaio relations are added, also from 
empirical values. The model galaxies are housed within dark mat¬ 
ter halos that follow the mass-concentration relation with scatter, 
as found in N-body simulations. We compare model populations 
that make different assumptions for dark matter halo profiles, and 
explore the consequent relations between Vrot(lkpc) and Vmax. 

When the model galaxy population is housed within halos 
with universal, steep, NFW density profiles, none of the galax¬ 
ies have low values of Vrot(l kpc)/Vmax, as found in some ob- 
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Figure 7. The scatter in the ratio Vrot(lkpc)A^max plotted as a function of 
Vmax for the SIDM model, shown for 1, 2 and 3 kpccores as pink, violet, 
and green triangles respectively, as well as for the observed galaxies (black 
circles). Error bars show the standard error: bins contain 11 (18) observed 
(model) galaxies. 


served galaxies, and the model population has too little scatter in 

Vrot(lkpc)/V max- 

The model population with the DC 14 density profile, whose 
profile is determined by the ratio M*/Mhaio, has galaxies with 
low values of Vrot(lkpc)/Vmax and a large scatter in the ratio 
Vrot(lkpc)/Vmax, both in reasonable agreement with observations. 
The model galaxies with low Vrot(lkpc)/Vmax result from halo ex¬ 
pansion, which is included in the DC 14 density profile. Two effects 
are responsible for the increased scatter compared with the NEW 
model: firstly, the increased fraction of baryonic mass in the inner 
regions when profiles are cored, which heightens the effect of the 
scatter in disc scale lengths. Secondly, the dependance of the DC 14 
density profile on M*/Mhaio means that halos with the same mass 
can have different core sizes. 

We note that the model that best fit observations included both 
halo expansion from non-adiabatic effects of outflows, as well as 
adiabatic contraction, both dependent on the ratio of M*/Mhaio. 
Both baryonic effects are required to explain the existence of galax¬ 
ies below and above the Vrot(lkpc)Vmax relation that is predicted 
by N-body simulations. 

We also looked at a specific case of a SIDM model, from 
which we can draw some general conclusions for models in which 
core size depends on Mhaio- The increased fraction of baryonic 
mass in the inner regions is also present in models where core 
size is dependent on halo mass. This means that the scatter in scale 
lengths results in such models having significantly more variation 
in rotation curve shapes than the NEW models. The cores of such 
models also result in low values of Vrot(lkpc)/Vmax in some galax¬ 
ies, in line with observations. Scatter in Vrot(lkpc)/Vmax depends 
on the core size, and can be as large as the observations for core 
sizes of 3kpc, using the definition of core size of eq. 3. However, 
core size may also have a mass dependance in SIDM models { |Vo-| 
Igelsberger et al.||2012| |, which is not accounted for in our simple 
parameterization. 

Eurther, it is likely that a scatter in core size versus Mhaio will 
exist in SIDM halos, with a dependence on merger histories. In¬ 
deed, [Rocha et al.| ( |2013| ) find scatter in core sizes at given halo 
mass. Euture studies will be able to quantify such scatter, and its 


mass dependence, which may result in a better match between the 
SIDM model and observations for smaller core sizes. It is certainly 
also possible (as recently shown in [Ery et al.|2015j l that baryonic 
physics would affect the dark matter density profiles in a SIDM 
universe, in a mass dependent manner. 

Our study suggests that baryonic physics are able to affect 
density profiles in dark matter models, to a degree that can account 
for the observed variation of rotation curve shapes. Such observa¬ 
tions may also provide constraints on alternative dark matter mod¬ 
els, such as SIDM, in terms of the scatter and size of cores at a 
given halo mass. 
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